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ABSTRACT 



Sr, Nd and Pb isotopic analyses of 477 samples representing 30 islands or 
island groups, 3 seamounts or seamount chains, 2 oceanic ridges and 1 oceanic 
plateau [for a total of 36 geographic features] are compiled to form a 
comprehensive oceanic island basalt [OIB] data set. These samples are 
supplemented by 90 selected mid-ocean ridge basalt [MORB] samples to give 
adequate representation to MORB as an oceanic basalt end-member. This 
comprehensive data set is used to infer information about the Earth's mantle. 
Principal component analysis of the OIB+MORB data set shows that the first 
three principal components account for 97.5% of the variance of the data. Thus, 
only four mantle end-member components [EMI, EMII, HIMU and DMM| are 
required to completely encompass the range of known isotopic values. Each 
sample is expressed in terms of percentages of the four mantle components, 
assuming linear mixing. There is significant correlation between location and 
isotopic signature within geographic features, but not between them, so 
discrimination analysis of the viability of separating the oceanic islands into those 
lying inside and outside Hart's (1984, 1988) DUPAL belt is performed on the 
feature level and yields positive results. 

A "continuous layer model" is applied to the mantle component percentage 
data to solve for the spherical harmonic coefficients using approximation 
methods. Only the degrees 0-5 coefficients can be solved for since there are only 
36 features. The EMI and HIMU percentage data sets must be filtered to avoid 
aliasing. Due to the nature of the data, the coefficients must be solved for using 
singular value decomposition [SVD], versus the least squares method. The F-test 
provides an objective way to estimate the number of singular values to retain 
when solving with SVD. With respect to the behavior of geophysics control data 
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sets, only the degree 2 spherical harmonic coefficients for the mantle components 
can be estimated with a reasonable level of confidence with this method. 

Applying a "delta-function model" removes the problem of aliasing and 
simplifies the spherical harmonic coefficient solutions from integration on the 
globe to summation over the geographic features due to the properties of delta- 
functions. With respect to the behavior of geophysics control data sets, at least 
the degree 2 spherical harmonic coefficients for the mantle components can be 
estimated with confidence, if not the degrees 3 and 4 as well. Delta-function 
model solutions are, to some extent, controlled by the nonuniform feature 
distribution, while the continuous layer model solutions are not. 

The mantle component amplitude spectra, for both models, show power at 
all degrees, with no one degree dominating. The DUPAL components [EMI, 
EMU and HIMU], for both models, correlate well with the degree 2 geoid, 
indicating a deep origin for the components since the degrees 2-3 geoid is 
inferred to result from topography at the core-mantle boundary. The DUPAL 
and DMM components, for both models, correlate well [and negatively) at degree 
3 with the velocity anomalies of the Clayton-Comer seismic tomography model 
in the 2500-2900 km depth range [immediately above the core-mantle boundary). 
The EMU component correlates well [and positively] at degree 5 with the 
velocity anomalies of the Clayton-Comer model in the 700-1200 km depth range, 
indicating a subduction related origin. Similar positive correlations for the geoid 
in the upper lower mantle indicate that subducted slabs extend beyond the 670 
km seismic discontinuity and support a whole-mantle convection model. 

Thesis Supervisor: Dr. Stanley R. Hart 

Title: Senior Scientist, Woods Hole Oceanographic Institution 
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CHAPTER 1 
INTRODUCTION 

PREVIOUS WORK 

That the Earth's mantle is heterogeneous is no longer a subject of 
controversy among geochemists, but the composition, the location and the 
geometry of these heterogeneities is very much in question. Direct sampling is 
not an option for studying the chemistry of most of the mantle, so products of 
indirect sampling, such as oceanic island basalts [OIB's] and mid-ocean ridge 
basalts [MORB], are invaluable for revealing the nature of the inaccessible 
mantle. Though the OIB's may be contaminated by interactions with the 
lithosphere or may sample large vertical sections of the mantle, they still retain 
the signature of their original source. 

Using various statistical methods and models, previous workers have 
defined what they believe to be the number of mantle component end-members 
required to represent the variation in the oceanic mantle data fOIB+MORB). 
Early on, Zindler el al. (1982) used factor analysis to evaluate the oceanic data in 
five dimensions. Their analysis indicated that the oceanic data define a plane [ the 
"mantle plane"], described by the mixing of three chemically independent 
components, two undifferentiated or slightly enriched mantle components and 
one MORB-type or depleted mantle component. 

Other workers have chosen five groups or components to represent the 
data. Using a series of two-dimensional isotopic plots. White (1985) divided the 
oceanic data into five distinct basalt groups [MORB, St. Helena, Kerguelen, 
Society, and Hawaii]. He concedes that the five groups may be end-members 
which mix to form intermediate compositions, but he believes that each group 
either represents a distinct, internally homogeneous reservoir or that each group 



is composed of a number of isotopically similar reservoirs. Likewise, Li el al. 

(1991) proposed fives extremes, using non-linear mapping: Atlantic MORB 
[DMM], St. Helena [HIMU], Walvis [EMI], Samoa [EMII] and D 5 [EMIII|. Non- 
linear mapping approximately preserves the geometric structure of the data by 
maintaining interpoint distances. Four of the five extremes of Li el al. (1991) 
are based solidly on samples trends from islands, but the D 5 extreme is based 
only on that one sample. More data is needed to substantiate their fifth extreme. 

By far the majority of analyses indicate the existence of four end-member 
components for the oceanic mantle data. Using two-dimensional plots, Zindler 
and Hart (1986) defined the following four end-member components: depleted 
MORB mantle [DMM], high U/Pb mantle [HIMU], and two enriched mantle 
components [EMI and EMII], with possibly two other components prevalent 
mantle composition [PREMA] and bulk silicate Earth [BSE]. Eigenvector 
analyses by Allegre el al. (1987) agree with the four component model of 
Zindler and Hart (1986). The four extremes of Allegre et al. (1987) are 
[correspond to]: extreme MORB [DMM]; St. Helena, Tubuai and MangaV islands 
[HIMU]; Kerguelen, Gough, Tristan da Cunha and Raratonga islands [ EM I ] ; and 
Sao Miguel and Atui islands [EMII]. Hart (1988), using an augmented data set 
and two-dimensional plots, concluded that the four end-members proposed by 
Zindler and Hart (1986) are valid representations of the extremes of the oceanic 
data. He resolves White's (1985) groupings into his own four component system 
as follows [White = Hart]: MORB = DMM, Society = EMII, St. Helena = HIMU, 
Hawaii = EMI, with the suggestion that White’s fifth group, Kerguelen, is a 
mixture of EMI and EMIL In addition, Li et al. (1991) also noted a tetrahedral 
structure to the data, when using factor analysis with varimax rotation, with the 
following four extremes: Atlantic MORB [DMM], Mangaia [HIMU], Samoa 
[EMII] and Walvis[EMI], 



One scenario for the genesis of the three unusual mantle components is put 
forth by Hart (1988). He proposes that HIMU, enriched in U, is probably 
generated by intra-mantle metasomatism, that EMI corresponds to a slightly 
modified bulk-earth compositon and that EMII can be explained by the recycling 
of sediments during subduction. The proposed formation mechanisms in no way 
limit the geometry of the mantle needed to generate the heterogeneities and, as 
such, a wide variety of models have been proposed. A whole mantle convection 
model might portray the enriched mantle components as blobs floating around in 
a depleted mantle matrix (Zindler and Hart, 1986) or perhaps as an accumulated 
layer of subducted oceanic crust and sediment at the core-mantle boundary that 
reaches the surface in mantle plumes (Hofmann and White, 1982). A two-layer 
convection model might rely on a depleted upper mantle feeding the mid-ocean 
ridges and an enriched lower mantle feeding oceanic islands via mantle plumes 
(Dupre and Allegre, 1983) or require a depleted upper mantle, a primitive lower 
mantle and an accumulated layer of subducted oceanic crust and sediment at the 
670 km discontinuity that supplies the enriched components via mantle plumes 
(White, 1985; Allegre and Turcotte, 1985). Anderson (1985) even proposes a 
three-layer convective model with the geochemical contrasts occurring only in 
the upper mantle with a depleted lower part that supplies the mid-ocean ridges 
and an enriched upper part from subduction of oceanic crust and sediment. 

A deep origin for the enriched components is indicated by Hart's (1984) 
large-scale isotopic anomaly, the DUPAL anomaly, characterized by the 
concentration of the enriched mantle components in a band from 2° S to 60° S. 
Qualitatively, countours of the anomaly criteria [A7/4, A8/4 and ASr (Hart, 

1984)] correspond to long-wavelength (and thus deep] geophysical quantities 
(Hart, 1988). Other researchers oppose this deep origin interpretation, citing the 
nonuniform distribution of hotspots as the reason for the pattern (White, 1985) 



or arguing that the DUPAL compositions occur in scattered locations and do not 
cover a coherent geographic area (Allegre et al., 1987). 

The purpose of this thesis is three-fold: (1) to address once again the issue 
of the number of mantle end-member components needed to represent the 
oceanic mantle data, (2) to statistically test the viability of the DUPAL distinction 
as a means of characterizing the OIB data and (3) to try to pinpoint the source 
depth of the enriched mantle components by expanding their relative abundances 
in spherical harmonics and comparing their expansions to those of known 
geophysical quantities. 



DATA 

The majority of this study focuses on Sr, Nd and Pb isotopic analyses of 
volcanic rocks from oceanic islands, seamounts, ridges, and plateaus. All of 
these geographic features overlie oceanic crust, with the exception of Nunivak 
Island on the Alaskan Continental Shelf, and none of them is directly associated 
with seafloor spreading, with the exception of Iceland, which has a mixture of 
mid-ocean ridge and hotspot influences. Essentially, the data set is that compiled 
by Zindler et al. (1982) and later augmented by Hart (1988), with some 
additional recent analyses (Appendix). Samples in the data set are mainly basalt, 

with some gabbros and trachybasalts; trachytes and other silica-rich rocks 
relative to basalt [roughly SiC >2 > 50% ] are excluded. The majority of the 

samples are of Cenozoic age, with the exception of the Walvis Ridge, Rio Grande 
Rise and New England Seamounts samples, with ages up to 100 Ma. If a choice 
is given, analyses of leached samples are preferred over analyses of unleached 
samples. In addition, only single samples for which there are Sr, Nd and Pb 
analyses are included. For consistency, Sr data is adjusted to 0.70800 [E&A 
standardl or 0.71022 (NBS SRM 987 standard] and Nd data is adjusted to 



0.512640 [BCR- 1 standard] or 0.51 1862 [La Jolla standard) or 0.51 1296 [Spex 
standard]. 

In this data set, referred to as the OIB data set, there are 477 samples 
representing 30 islands or island groups, 3 seamounts or seamount chains, 2 
aseismic oceanic ridges and 1 oceanic plateau (Figure 1.1 and Table 1.1). The 
isotopic means and standard deviations for the OIB data are listed in Table 1.2. 

Since MORB is considered to be one of the mantle component end- 
members (Zindler el al., 1982; White, 1985; Zindler and Hart, 1986; Hart, 

1988), any attempt to choose end-members should include MORB data. For this 
reason, a second data set is created using the OIB data and a selection of 90 
MORB samples (Appendix), the OIB+MORB data set (Table 1.3). The criteria 
for choosing OIB samples applies to the MORB samples as well. Isotopic means 
and standard deviations for the OIB+MORB data are listed in Table 1.2. 

ORGANIZATION 

The main thust of this work is to characterize the OIB data and to search 
for possible correlations between the geochemical signatures of OIB's and 
geophysical quantities, such as the geoid and seismic tomography, that might help 
pinpoint the depth[s] of the OIB reservoirfs]. 

Chapter 2 explores the nature of the OIB isotope data. With the help of 
principal component analysis, the data is expressed in terms of percentages of 
four mantle component end-members. Spatial correlation testing reveals the 
relationship between geographic distance from island to island and feature to 
feature and the "isotopic distance" between samples. Discrimination analysis, 
both nearest-neighbor and graphical, is used to test the viability of separating the 
oceanic islands into two groups, inside and outside the DUPAL belt. 



Chapter 3 applies a "continuous layer model" to the mantle component 
data, as an assumed geometry for the OIB reservoir, in order to solve for the 
spherical harmonic coefficients. The problem of aliasing is addressed with the 
relationship of variation in mantle components to distance between features. 
Approximation methods are used to solve for the coefficients. Geophysical data 
sets are constructed, using GEM-L2 geoid coefficients, to serve as controls 
against which to judge the success of the approximation methods. 

Chapter 4 applies a "delta-function model" to the mantle component data 
to provide a mathematically more robust solution for the spherical harmonic 
coefficients. The delta-function approximation removes the problem of aliasing, 
but generates a solution dependent upon feature location. The same geophysical 
data sets are used again to judge the success of the delta-function approximation. 

Chapter 5 compares the mantle component spherical harmonic solutions 
for the two models in terms of their amplitude spectra, how well they correlate 
with the geoid, how they are affected by the nonuniform feature distribution and 
how well they correlate with the Clayton-Comer seismic tomography model. 

The implications of these results and recommendations for further research are 
discussed. 



Table 1.1. Geographic features represented in the 01B data set, with their 
components, number of samples (in braces] and references indicated. 









Feature 


Components 


References * 


Ascension [5] 




7,34,35 


Amsterdam/St. Paul [11] 




Amsterdam [5] 


38 




St. Paul [6] 


38 


Azores [6] 




Faial [1] 


22 




Sao Miguel [5] 


L8 


Balleny [3] 




19 


Cameroon Line [18] 




Bioko [5] 


17,18 




Pagalu [ 1 ] 


18 




Principe [3] 


18 




Sao Tome [9] 


17,18 


Cape Verde Islands [41] 




Fogo [6] 


14 




Maio [9] 


8,14 




Sao Antao [10] 


8,14 




Sao Tiago [13] 


14 




Sao Vincente [3] 


14 


Christmas [13] 




19 


Cocos [3] 




3 


Comores Archipelago 1 14] 




38 


Cook-Austral Islands [26] 




Aitutaki [4] 


1,21 




Atui [6] 


1,21 




Mangia [5] 


1,21 




Mauke [3] 


1.21 




Raratonga [8] 


1,21,23 


Crozet Islands [9] 




38 


Fernando de Noronha [16] 




1,13 


Galapagos Islands [11] 




39 


Gough [2] 




1 


Hawaiian Islands [73] 




Hawaii [14] 


28,32 




Kahoolawe [13] 


37 




Kauai [2] 


28 




Lanai [4] 


37 



Table 1.1. Continued. 



Feature 


Components 


References 


Hawaiian Islands [73] 




Loihi [15] 


27 




Maui [3] 


28 




Molokai [5] 


28 




Oahu [17] 


29 


Iceland [7] 




20 


Juan Fernandez Islands [4] 




15 


Kerguelen Plateau [41] 




Heard Island [9] 


2,30 




Kerguelen Island [20] 


12,30,38 




Kerguelen Plateau [12] 


26,36 


Louisville Seamount Chain [4] 




6 


Marion/Prince Edward [4] 




19 


Marquesas Archipelago [11] 




10,11,33 


Mascareignes [8] 




Mauritius [1] 


1 




Reunion [7] 


38 


New England Seamounts [6] 




31 


Nunivak [2] 




25 


Pitcairn [19] 




41 


Ponape [ 1 ] 




19 


Sala Y Gomez [1] 




? 


Samoa Islands [34] 




Manu'a [4] 


42 




Savai’i [8] 


42 




Tutuila [9] 


23,42 




Upolu [13] 


23,42 


San Felix/San Ambrosio [5] 




San Felix [4] 


15 




San Ambrosio [1] 


15 


Shimada Seamount [ 1 ] 




16 


Society Ridge [9] 




Mehetia [2] 


9 




Moua Pihaa [ 1 ] 


9 




Tahaa [1] 


40 




Teahitia [4] 


9 




dredge [ 1 ] 


9 


St. Helena [31] 




7,1,4,7.22 



Table 1.1. Continued. 



Feature 


Components 


References 


Trinidade [1] 




1 


Tristan de Cunha [5] 




7,22 


Tubuai-Austral Islands [22] 


Marotiri [1] 


5 




Raevavae [ 1 ] 


1 




Rapa [3] 


5,23 




Rimatara [4] 


?,21 




Rurutu [4] 


21,23 




Tubuai [9] 


5 


Walvis Ridge [10] 




24 



Mn the reference column, a "?" indicates a sample with an unknown reference. 

Reference guide: [1] AUegre et al., 1987; [2] Barling and Goldstein, 1990; |3| 
Castillo et al., 1988; [4] Chaffey et al., 1989; [5] Chauvel et al., 1991; 1 6 J Cheng 
et al., 1988; [7] Cohen and O'Nions, 1982a; [8] Davies et al., 1989; |9] Devey et 
al., 1990; [10] Duncan et al., 1986; [1 1 j Dupuy et al., 1987; [12] Gautier et al., 
1990; [13] Gerlach et al., 1987; [14] Gerlach et al., 1988; [15] Gerlach el at., 
1986; [16] Graham, 1987; [17] Halliday etal., 1990; [18] Halliday et al., 1988; 
[19] Hart, 1988; [20] Hart, unpublished; [21] Nakamura and Tatsumoto, 1988; 

[22] Newsom et al., 1986; [23] Palacz and Saunders, 1986; [24] Richardson et al., 
1982; [25] Roden, 1982; [26] Salters, 1989; [27] Staudigel et al., 1984; [28] Stille 
et al., 1986; [29] Stille etal., 1983; [30] Storey et al., 1988; [31] Taras and Hart, 
1987; [32] Tatsumoto, 1978; [33] Vidal etal., 1984; [34] Weis, 3983; [35] Weis et 
al., 1987; [36] Weis et al., 1989; [37] West et al., 1987; [38] White, unpublished; 
[39] White and Hofmann, 1982; [40] White etal., 1989; |41] Woodhead and 
McColloch, 1989; [42] Wright and White, 1987. 



